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Abstract 

An important radiation field for photomeson neutrino production in blazars is shown to be the 

o i 

radiation held external to the jet. Assuming that protons are accelerated with the same power as 

^— > . 

electrons and injected with a —2 number spectrum, we predict that km 2 neutrino telescopes will 
\q ■ detect <^ 1 neutrinos per year from flat spectrum radio quasars (FSRQs) such as 3C 279. The 

, escaping high-energy neutron and photon beams transport inner jet energy far from the black-hole 

> 

fT) ' engine, and could power synchrotron X-ray jets and FR II hot spots and lobes. 
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EGRET observations of ~ 100 MeV - 5 GeV emission from over 60 FSRQs and BL 
Lac objects established that blazars are, with gamma-ray bursts, among the most powerful 
accelerators of relativistic particles in nature [|1|. The standard blazar model consists of a 
supermassive black hole ejecting twin jets of relativistic plasma, one of which is pointed 
towards us. Observations of bright, highly variable nonthermal radiation from blazars mean 
that regions of intense photon and nonthermal particle energy densities, both of which are 
needed for efficient photopion production |2|, Hf , are found in these sources. Blazars therefore 
represent a potential source of energetic 7r ± -decay neutrinos to be detected by operating and 
planned high-energy neutrino telescopes Q. Previous treatments have considered internal 
synchrotron photons and the direct disk radiation field || as targets for high-energy proton 
interactions in jets, and neutrino production in the cores of AGN (active galactic nuclei) ||. 
Here we show that photons from external quasi-isotropic radiation fields, which have earlier 
been proposed as target photons to be Compton-scattered by nonthermal electrons to 7-ray 
energies 0, also provide the most important photon source for photomeson production of 
^ 30 TeV neutrinos in FSRQs. The inclusion of this effect increases neutrino detection 
rates by more than an order-of-magnitude. Moreover, the neutrinos are formed at energies 
^3 x 10 13 eV rather than at ^ 10 17 eV [§], which improves prospects for detection. 

Strong optical emission lines || from the illumination of broad-line region (BLR) clouds 
reveal bright accretion-disk and scattered disk radiation in the inner regions of FSRQs. 
BL Lac objects have weak emission lines, so the dominant soft photon source is thought 
to be the internal synchrotron emission in BL Lac objects. In our analysis of photomeson 
production in FSRQs, we assume that the quasi-isotropic scattered external radiation field 
dominates the direct accretion-disk field. The photomeson neutrino spectrum can be calcu- 
lated once the mean magnetic field and comoving spectral energy density are determined. 
The measured variability time scale t var and synchrotron (L s ) and Compton (Lc) luminosi- 
ties determined from the spectral energy distribution of the bright, well-studied blazar 3C 
279 |Hj] implies physical parameters of its jet emission region. Its redshift z = 0.538, which 
implies a luminosity distance g?l — 1-05 x 10 28 cm for an Q m = 0.3, = 0.7 cosmology 
with a Hubble constant of 65 km s _1 Mpc -1 . A crucial unknown is the Doppler factor 
5 = [r(l — /?rcos^)] _1 , where T is the bulk Lorentz factor of the relativistic plasma blob, 
and 9 is the angle between the jet and observer directions. 

The comoving synchrotron photon energy density is given by u' s = L s / '(27rr 2 c<5 4 ) , where 
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Tf, is the comoving radius of the blob, here assumed spherical, and primes denote comov- 
ing quantities. The spectral energy density u s (e') = m e c 2 e'n' s (e') (units of ergs s^V -1 ) of 
photons with dimensionless comoving energy e' = hv'/m e c 2 is found through 

, ~ l<tjfs{e) „ 2^(1 + zff s {e) 

eUs[e) ~ rlc8* ~ <?tl ar 8* ' K } 

in terms of the vF v flux / s (e) of the synchrotron component. Here e = Se'/(1 + z), and we 
relate t var and through the expression r b = ct var 5/(l + z). 

For calculations of n' s (e') from 3C 279, we approximate the flux density F(e) oc e~ a 



observed during the flare of 1996 |Tl|] in the form of a continuous broken power-law function, 
with indices a.\ = 0.5, c*2 — 1.45, and = 0.6 at frequencies z/ < v < u\ = v p k, with 
v p k — 10 13 ^i3 Hz, v \ < v < v 2 = 10 16 Hz, and v > i>2, respectively. The vF v synchrotron 
radiation flux f s (e) reaches a maximum value f s {t P k) — 1-7 x 10 -10 ergcm -2 s _1 at v = v v ^ 
or e = e p k- The flaring and three- week average flux of the Compton 7-ray component peaks 
at ~ 500 MeV and is ~ 20 and ~ 5-10 times larger than f s {e p k) at the corresponding times 
|n|. Using t var (d) ~ 1 (in days) for 3C 279 as observed by EGRET [[Tl]] during the three 



week average, one finds from equation (JI|) with <5i = 5/10 that 

< ^ e' pk u' s {e' pk ) ~ 0.4 [t var (d)]- 2 ^o 6 erg cm" 3 . (2) 

If the 7-rays from FSRQs are due to Compton-scattered radiation from external photon 
fields, then L EC /L S = u' ext /uB 0, where L EC = L c — L ssc is the measured power from 
Compton-scattered external photon fields, and Lssc is the synchrotron-self- Compton (SSC) 
power. Consequently u' ext = ub(L ec /L s ) = au' s = aL s / ^(27rr 2 c5 4 ), where a = L EC / L S sc- 
The energy of the external photons in the comoving frame is e' ext = Se ex t- Models taking 
into account SSC and external Compton (EC) components show that a complete spectral 



fit requires synchrotron, SSC and EC components [T2], with a ~ 0.1-1 for BL Lac objects 
and a ~ 1-10 for FSRQs. During the 1996 flare of 3C 279, a ~ 10, so that the energy 
density of external soft photons in the comoving blob frame is ~ 10 x greater than the 
synchrotron energy density. Moreover, the energy of these photons increases with 5, lowering 
the photomeson production threshold. 

Lower limits to 5 are defined by the condition that the emitting region be transpar- 



ent to 77 pair-production attenuation inside the blob |]13 |, which requires that T 77 (e') = 
2<TT^ / s (2/e')rfc/(3e / ) < 1. The measured > 10 16 Hz flux during the 1996 flare then implies 
that 5 > 5[£ M (GeV) ai 7[t mr (d)] a19 , where E ph (GeV) is the photon energy in GeV. 



The magnetic field B in the blob can be determined by introducing the equipartition 
parameter rj = u' el (l + k pe )/u-B for the ratio of relativistic electron to magnetic energy 
densities in the jet, with the factor k pe correcting for the presence of nonthermal hadrons. 
We assume k pe = 1 in the estimate of B. The measured synchrotron flux density F s oc e~ a 
with a ~ 0.5 at v < 10 13 Hz gives the equipartition magnetic field 

5 (Gauss) 2 130 ^[(l + ^)»n(V^(l + z)^ 

where vq/vx = 10 3 and /_ 10 = / s (e pfe )/(10- 10 erg cm" 2 s^ 1 ) = 1.7 for 3C 279. 

An alternative estimate of B can be derived from the observed synchrotron and Comp- 
ton powers |TJ]]. Neglecting Klein-Nishina effects on Compton scattering, the ratio of the 



Compton and synchrotron powers Lc/L s = u' ph /uB, so that ub = (L s /Lc)u f s (l +u' ext /u' s ) 
(L s /Lc)u' s (l + a). Using equation (TO) in this expression for B gives 



2(1 + 2) WTTfl 

The observed flux of synchrotron radiation from 3C 279 with a = 0.5 at v <C 10 13 Hz 
implies that B(G) ~ 7 r/ - 2//7 [t var (d)]~ 6 / 7 <5 10 13//7 from equipartition arguments. The estimate 
of B from equation (§ gives B(G) = 4 t^^d)^ 3 , using L s = 3 x 10 47 ergs s -1 , L c /L s = 10, 
and a = 10. The external radiation energy density in the comoving frame, using B derived 
from equation (|^), is u' ext ~ 1.5 {Lec/ L s )rj~^ 7 \t vax {d)Y 12 / 7 eig cm -3 . Note the weaker 
dependence of u' ext on 5 compared to u' s . For the calculations we take r\ = 1 and L EC /L S = 10. 

Energy losses of relativistic protons (and neutrons) are calculated on the basis of standard 
expressions (e.g., Ref. fll5| ) for the cooling time of relativistic protons due to photopion 
production in jry collisions. If the ambient photons have spectral density n' h(e'), then 



27p 2 e /: 

where 7 P is the proton Lorentz factor, e r is the photon energy in the proton rest frame, cr(e r ) is 
the photopion production cross section, e t h ~ 150 MeV is the threshold energy for the parent 
photon in the proton rest frame, and K vl {e T ) is the inelasticity of the interaction. The latter 
increases from K\ pa 0.2 at energies not very far above the threshold, to K2 ~ 0.5 — 0.6 at 



larger values of e r where multi-pion production dominates |15|, 

A detailed recent study of this photohadronic process is given by Ref. ||16|| . To simplify 
calculations, we approximate cr(e r ) as a sum of 2 step-functions o"i(e r ) and cr 2 (e r ) for the total 



single-pion (p+7 — > p+7r° or n+7r + ) and multi-pion channels, respectively, with o\ = 380 fib 
for 200 MeV < e r < 500 MeV and a% — outside this region, whereas 02 = 120 /ib at 
e r > 500 MeV. The inelasticity is approximated as K pi = K\ and K pi = K 2 below and above 
500 MeV. Our calculations give good agreement with more detailed treatments of the time 
scales for photopion interactions of ultra-relativistic cosmic rays with the cosmic microwave 
background |T5|, [T7| . This approach also works well for a broad power-law distribution of field 
photons u' h (e') oc (e')~ a for different spectral indices a, and readily explains the significant 
increase in the mean inelasticity of incident protons (or neutrons) from (K P1 ) ~ 0.2 for steep 
photon spectra with a 7 > 1, to (K P1 ) — > 0.6 for hard spectra with a 7 < 1 ff6|. 

The spectra of secondary 7r 0,± -decay particles (V, 7, e) are calculated in the ^-function 
approximation, assuming that the probabilities for producing pions of different charges (n , 
7r + and 7r~) are equal for the multi-pion interaction channel. To correctly apply the 5- 
function approximation, one has to properly take into account the different inelasticities of 



the multi-pion and single-pion production channels (see Ref. [16] for comparison). 

We assume that the spectrum of the external UV radiation field arises from a 
Shakura/Sunyaev [jTBl optically-thick accretion disk model that is scattered by BLR clouds, 
though other sources of external photons could be reprocessed line emission or infrared ra- 
diation from a torus. This disk model has flux density F(e) oc e 1//3 exp(— e/e max ), where the 
maximum photon energy e max is determined by the innermost radius of the blackbody disk 
and properties of the central engine. We take m e c 2 e max = 35 eV . 

Fig. 1 shows observer-frame energy loss timescales of protons due to photopion production 
in a jet of 3C 279 calculated for Doppler-factors 5 = 7, 10, and 15. For comparison, we show 
also the photopion timescales if only the interactions with the synchrotron radiation were 
taken into account for the cases 5 = 7 and 10. Fig. 2 shows the energy distributions of 
relativistic protons N p (E) in the jet plasma blobs calculated assuming power-law injection 
of relativistic protons with number index a p = 2 on observed timescales At = 3 weeks. 
The total injection power of the protons L p = 2 x 10 48 <5~ 4 erg s -1 , which corresponds to 
the average measured 7-ray luminosity from 3C 379 during the three week observing period 
surrounding the 1996 flare |[L1]| . In calculations of iVp we take into account the photohadron 
interaction energy losses, as well as the escape losses of the protons in the Bohm diffusion 
limit. This limit ensures that the gyroradii rL(cm) = 3.1 x 10 6 (A7 P /Z)/B(G) of the highest 
energy ions do not exceed r& . 
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Importantly, we also take into account proton losses through the p7 — > nX channel. 
The neutrons may then escape the blob either before they decay or collide again with 
photons inside the blob. To demonstrate the significance of this channel, in Fig. 2 by the 
full dotted curve we show the proton distribution calculated for the case of 5 = 7, but 
neglecting the effect of escaping neutrons. Comparison with the solid curve indicates that 
~ 50% of the E > 10 14 eV proton energy is carried by relativistic neutrons from the blob, 
forming a highly collimated neutron beam with opening angle 9b ~ 1/T and a total power 
L (beam) _ Q.25L p ~ 5 x 10 47 <5 4 ergs 1 for the parameters assumed in Fig. 2. 

Fig. 3 shows the expected differential fluences of neutrinos produced by protons in Fig. 
2, integrated over the 3 week flaring period. Note that the mean neutrino energy fluxes 
at E > 100 TeV are at the level E F V (E) ~ E 2 $ v /At < 2 x 10" 11 erg cm" 2 s -1 . Given 



very similar emissivities in overall secondary v and (7 + e) production (e.g. Ref. |16|), the 
unabsorbed fluxes of multi-TeV 7-rays would therefore be at the same level, which makes 
only < 2% of the 7-ray energy flux detected by EGRET. A pair-photon cascade in the blob 
would increase this fraction to the level of < 10% contribution of the cascade radiation to 
the detected 7-ray flux, reprocessing most of the initial 7-ray energy into 7-radiation with 
E ^5 10 GeV. A small fraction of the energy will appear as synchrotron radiation in the 
X-ray to GeV 7-ray band. Because the cascade synchrotron power is emitted with a very 
hard spectrum over a wide energy range, it will not overproduce also the observed X-ray flux 
of 3C 279. These estimates show, however, that the nonthermal proton power cannot be 
more than 10 x the 7-ray power, or the cascade radiation would exceed the observed fluxes. 

For the fluences in Fig. 4, the total number N v of neutrinos that could be detected by a 
lkm 2 detector, using neutrino detection efficiencies given by Ref. [|J, is 0.45, 0.27 and 0.12 
for the cases of 5 = 7, 10 and 15, respectively. These numbers are, however, much less if the 
external field is neglected, giving instead the values N v ~ 0.055 (5 = 7) and 0.014 (5 = 10) 
for the fluences in Fig. 3. This would not leave a realistic prospect for the detection of at 
least 2-3 neutrinos, which is required for a positive detection of a source. BL Lac objects, 
which have weak BLRs and, therefore, a weak scattered external radiation field, should 
consequently be much weaker neutrino sources. 

For a 10% flaring duty factor, and considering the additional neutrino production during 
the quiescent phase, we expect that IceCube or other km 2 array may detect several neutrinos 
per year from 3C 279-type blazar jets with 5 ~ 5-15. Allowing an overall power of protons 
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larger than that for primary electrons improves the prospects for neutrino detection. A 
softer nonthermal proton spectrum would, however, reduce the neutrino emissivity. 

Acceleration of protons only to E p > 10 14 eV is required for efficient neutrino production 
through photomeson interactions, as can be seen from Fig.l. Thus we predict that a km 2 
array will detect high-energy neutrinos from FSRQs without requiring acceleration of protons 
to ultra-high (> 10 19 eV) energies. The maximum cosmic-ray energies from blazars are found 
to be EcR,max — 8 x 10 17 ZT 2 B (G)t mr (d) / (1 + z) eV p0| , with maximum neutrino energies 
a factor ~ 20 lower. Given additionally the intense photomeson interactions encountered by 
a proton at this optimistic limit in the jet environment (see Fig. 1), it is unlikely that the 
inner jets of blazars are significant sources of ultra-high energy cosmic rays. The neutrino 
production from FSRQs is therefore unaffected by the cosmic ray bound of Ref. (though 
note the more general upper limit derived in Ref. 0). 

About half of the photomeson interactions in the plasma jet will produce neutrons. Neu- 
trons with comoving Lorentz factors 7^ > 10 3 5iot var (d) will escape the blob in a collimated 
beam, and continue to interact efficiently with the external UV radiation field up to ~ 0.1 — 1 
pc scales, initiating pair-photon cascades far away from, but in the same direction with, the 
inner jet due to secondary 7-rays and electrons. FSRQ jets will thus form intense neu- 
tron and photon beams that transport energy far from the central engine of the AGN . 
Cosmic-ray neutrons with energies of 10 17 Eij eV will decay at distances ^ En kpc from the 
jet core. The deposition of the energy of the beamed hadrons and multi-TeV 7-rays (through 
77-absorption at > 10 kpc scales) can reasonably be at the mean rate ~ 10 43 ergs _1 , which 
could lead to enhanced X-ray synchrotron emission along the large-scale jet, as seen by the 
Chandra X-ray Observatory from sources such as Pictor A p3| . Depending on the stability of 



charged particle transport, nonthermal particles could be deposited at ^ 100 kpc to supply 
power for hot spots and lobes in FR II galaxies, though in situ particle acceleration could 
also take place at these sites. The predicted weaker neutron and neutrino production in BL 
Lac objects and their parent FR I galaxies might account for the morphological differences 
of these radio galaxy types. This picture will be tested by sensitive high-energy neutrino 
detectors, ground-based > 50 GeV gamma-ray detectors, and GLAST ||, |2~3}| . 
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FIG. 1: Photomeson interaction energy loss time scale of protons, calculated for spectral fluxes 
observed from 3C 279 (see text) and t var = 1 d, assuming 3 different Doppler-factors for the jet: 
5 = 7 (solid curve), 5 = 10 (dashed curve), and 5 = 15 (dot-dashed curve). The dotted and 
triple-dot-dashed curves are calculated for S = 7 and 5 = 10, respectively, when pj interactions 
with the synchrotron radiation field alone are considered. 
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FIG. 2: The energy distribution of relativistic protons, calculated for the same jet Lorentz-factors 
5 = 7, 10, and 15 and assumptions for the external and internal radiation fields as made in Fig. 1, 
assuming a power-law injection of relativistic protons with number index a p = 2 during At = 3 
weeks with a power L p = 2 x 10 48 5~ 4 ergs s _1 . The full dots show the proton distribution calculated 
without the effect of neutrons escaping the blob. 
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FIG. 3: The fluences of neutrinos expected due to photomeson interactions of protons shown in 
Fig. 2 with the external (UV) and internal (synchrotron) photons in the jet of 3C 279 for Dopper 
factors 5 = 7 (solid curve) 5 = 10 (dashed curve), and 8 = 15 (dot-dashed curve). The dotted 
and triple-dot-dashed curves show the neutrino fluences if the external photon field in the jet is 
neglected. 



12 



